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ABSTRACT 

We develop a method to analyze the effect of an asymmetric supernova on hierarchical 
multiple star systems and we present analytical formulas to calculate orbital param- 
eters for surviving binaries or hierarchical triples and runaway velocities for their 
dissociating equivalents. The effect of an asymmetric supernova on the orbital param- 
eters of a binary system has been studied to great extent (e.g. iHillsl 119831 : iKalogeral 
Il996t iTauris &: Takenj Il998l) , but this effect on higher multiplicity hierarchical sys- 
tems has not been explored before. With our method, the supernova effect can be 
computed by reducing the hierarchical multiple to an effective binary by means of 
recursively replacing the inner binary by an effective star at the center of mass of that 
binary. We apply our method to a hierarchical triple system si milar to the progenitor 
of PSR J1903+0327 suggested bv lPortegies Zwart et all (|201l[ ). We confirm their ear- 
lier finding that PSR J1903-I-0327 could have evolved from a hierarchical triple that 
became unstable and ejected the secondary star of the inner binary. Furthermore, if 
such as system did evolve via this mechanism the most probable configuration would 
be a small supernova kick velocity, an inner binary with a large semi-major axis, and 
the fraction of mass accreted onto the neutron star to the mass lost by the secondary 
would most likely be between 0.35 and 0.5 

Key words: stars: supernovae: general - binaries: general - methods: analytical - 
methods: numerical - pulsars: individual: J1903-I-0327. 



; 1 INTRODUCTION 

Asymmetric supernovae (SNe) in binary and hierarchical 
multiple star systems form a crucial phase in the forma- 
tion of stellar systems containing a compact stellar remnant 
- neutron star or black hole. In previous studies of SNe in 
binaries, two effects of the SN have been considered: (1) 
sudden mass loss and (2) a random kick velocity imparted 
on the compact remnant of the star undergoing the SN. The 
combined effect which changes the orbital parameters causes 
the binary to dissociate in the majority of the cases. 

The study of binaries surviving a supernova (S N) explo- 
sion o f one of its components was first performed bv lBlaauwl 
(|l96ll ') anddoersma (1961), assuming a symmetric SN (i.e. 
only mass loss). The necessity of asymmetry in the SN, re- 
sultin g in the kick velocity, was first suggested bv lShklovskiil 
The stat i stical study on pulsar scale heights by 
iGunn fc Ostrike"il (|l97Cll ) firmly supported the asymmetric 
SN model and to date the adding of the kick velocity to the 
newly born neutron star (or black hole) is a commonly ex- 

* The first and second authors contributed equally to the work. 



cepted mechanism (|van den Heuvel &: van Paradiis|[T997l ). 
Both the type of explosion mechanism and whether the ex- 
ploding star is in a binary syste m are found to influence th e 
effect of the kick velocity (see e.g. lPodsiadlowski et al. Il2004l ). 
but the exact physical process underlying the production of 
kicks remains unclear. The analysis of the effect of asym- 
metric supernovae on binaries has been sufficient to explain 
most of the observed post-SN stellar systems, and little to 
no effort has gone into studying the effect on hierarchical 
multiple star systems. 

Millisecond pulsar ( MSP) J1903-I-0327 (spin period ~ 
2.15 ms), first observed byfCh ampion et all (|2008l ) and later, 
in more detail, by iFreire et al . (2011), is part of what may 
be the first observed MSP binary to have evolved from 
a hierarchical triple progenitor. MSP J1903+0327 is or- 
bited by a main sequence star in a wide (orbital period 
~ 95.2 days) and eccentric (eccentricity e ~ 0.44) orbit. 
Based on these observables it seems impossible that this 
binary (hereafter J1903-f0327) for med via the traditiona l 
mechanism in a binary prog enitor (|Champion et al]|2008l ). 
iPortegies Zwart et all l|201ll ) proposed that the progenitor 
system was a binary accompanied by a third and least mas- 
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sive main-sequence star in a wider orbit about this binary. 
During the low-mass X-ray binary (LMXB) phase of the in- 
ner binary, the orbit of the LMXB expanded due to mass 
transfer from the evolving inner companion (donor) star to 
the neutron star, which was formed in the SN. This eventu- 
ally caused the triple to become dynamically unstable and to 
eject the inner companion resulting in the observed system 
J1903-h0327. 

J1903-f0327 is not a unique case, however: there is 
a significant number of sys tems like the progen i tor o f 
J1903-I-0327 as suggested in IPortegies Zwart et al.l (|201ll ') 
and similar hierarchical stellar systems of higher multiplic- 
ity. The Multiple Star Catalog lists 602 triples, 93 quadru- 
ples, 22 quintuples, 9 sextuples and 2 septuples l|Tokovininl 
[1993) of which 90 systems contain at least one star with a 
mass M ^ 10 Mq. Each of these multiples will eventually 
experience a core-collapse SN of the most massive star. After 
the SN these systems are either fully dissociated, dissociate 
into lower multiplicity multiple star systems, or survive the 
SN. 

We begin the study of the effect of an asymmetric SN on 
hierarchical multiple star systems by first readdressing the 
SN effect on a binary and subsequently treating the effect 
in a hierarchical triple. We show that a hierarchical triple 
can effectively be regarded as a binary system comprised 
of the center of mass of the inner binary and the tertiary 
star. The effect of a SN on a hierarchical triple system, now 
reduced to an effective binary, can be calculated using the 
prescription for a SN in binary. We ultimately generalize this 
effective binary method to hierarchical multiple star systems 
of arbitrary multiplicity. In the second part of the paper we 
perform Monte Carlo simulations of a hierarchical triple star 
sy stem similar to the pr o genito r of J1903+0327 suggested 
in IPortegies Zwart et al.l l|201ll ) to determine the (stable) 
survival rates, and evaluate whether such a formation route 
is plausible. 



2 CALCULATION OF POST-SN 
PARAMETERS 

2.1 Binary systems 

We consider a binary system of stars with mass, position 
and velocity for the primary and secondary star, given by 
(mi,o,ri ,vi,o) and (m2,r2,V2,o) respectiveljQ, in which the 
primary undergoes a SN. The binary system is uniquely de- 
scribed by the semi-major axis, ao, eccentricity, eo, and true 
anomaly, 60 ■ The separation distance is ro. We assume that 
the SN is instantaneous, meaning an instantaneous removal 
of mass of the primary, no SN-shell impact on the compan- 
ion (secondary) star, and the orbital motion during this mass 
loss phase is neglected, i.e. r = ro and V2 = V2,o. 

After the SN the orbital parameters have changed to: 
semi-major axis, a, eccentricity, e, and true anomaly, 9. For 
a general Kepler orbit of two objects with masses mi and 



^ Tlie contingent sufBx 1,2, etc. indicates which star we are con- 
sidering (e.g. 1 for the primary). The contingent suffix denotes 
the pre-SN state and when it is absent, it either refers to the 
post-SN state or the absence indicates that there is no difference 
in the pre- and post-SN states of that parameter. 



m2 respectively, a relative velocity, v, semi-major axis, a, 
and separation distance, r, the orbital energy conservation 
equation is 



G(mi 



- m2. 



(1) 



where G is Newton's gravitational constant. The specific 
relative angular momentum h is related to the orbital pa- 
rameters as follows 



Ir X vl' 

G(mi -|- m2)a(l 



(2) 
(3) 



where the first equality holds for all Kepler orbits and the 
second only applies to bound orbi ts. Fo r thor o ugh stud- 
ies on SNe in a binary syst e m see iHillj (|l983l ). iKalogeral 
(|l996t ). and lTauris fc Takend (|l998t ): the latter authors also 
take into account the shel l impact on the companion s tar us - 
ing a method proposed bv lWheeler. Lecar. fc McKed (|l975h . 
Following the mentioned works as guides for our calcula- 
tions on the binary system we use a total pre-SN mass of 
Mo = mi,o + rn2. Without loss of generality, we choose a 
coordinate system in which at f = the orbit lies in the 
xy-plane, the center of mass of the binary (cm) is at the 
origin, the y-axis is the line connecting the primary and the 
secondary (the cm coordinate system; see Figure [ij, and we 
choose a reference frame in which at t = the cm is at rest 
(the cm reference frame). 

Before the SN the separation distance between the stars 



r = ri - r2 = 0, - 



Using the following notation 



ao(l - Eq) 
1 -f- eo cos ^0 



(4) 



ao y'l — Bq cos 70 cos 60 + ao sin 70 sin 60 , 
— ao yl — ei cos 70 sin 60 + ao sin 70 cos 60 , 

X 



VOx 



VOy 



vo 



Vo 



y 



in which 70 is the pre-SN eccentric anomaly defined by r = 
ao(l — eo cos 70), the velocity of the primary relative to the 
secondary is 



Vo = Vl,o — V2 = {vox,vov,0)- 



(5) 



After the SN the primary has lost a part of its mass. Am, 
and has obtained a velocity kick Vk in a random direction, 
which makes an angle </> with the pre-SN relative velocity 
Vo. The velocity of the primary relative to the secondary, 
after the SN, is 



V = Vo +Vk = {vox +Vkx,VOy +Vky,Vkz), 



(6) 



the mass of the primary is mi — mi.o — Am and the total 
binary mass is M = Mo — Am. Applying these relations 
and equations ((l} and ([2]) to the binary system, we obtain 
equations relating the post-SN semi-major axis, a, and ec- 
centricity, e, to both the pre- and post-SN orbital parameters 
and velocities. Using «c,o = vo\r^ao = (GMo/ao)^^^ as the 
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a. The cm coordinate system in the cm reference frame for a 
binary system before the SN (at t = 0). 
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b. The cm coordinate system in the cm reference frame for a 
binary system after the SN. 

Figure 1. Schematic representation of a binary system in the 
pre- and post-SN phase. The solid blue circles denote the primary 
and secondary star; the solid red cirle denotes the cm. The solid 
arrows denote the velocities the stars or cm have at that phase; 
the dashed arrows denote the velocity the SN imposes on the 
stars or cm which will change its velocity in the next phase, a. 
In the pre-SN phase the coordinate system is centered on the cm 
being at rest. b. In the post-SN phase the coordinate system is 
no longer centered on the cm - the cm has been translated in 
the y-direction, towards the secondary, and has gained a velocity 
Vsys- In both cases the inner binary orbital plane lies in the xy- 
plane and the y-axis is the line connecting the primary and the 
secondary. 



pre- SN relative velocity (|Hillslll983l ). we obtain 



a 
ao 



2ao Am 
~ Mo 



„ "0 l^fc 

2 cos < 

Vc,0 Vc,0 



(7) 



eo) 



(Mo - Am) 2 



2ao Am 
~ Mo" 



-2 cos ( 

Vc,0 Vc,0 

a§(l-eg)2 



(8a) 



{-"Ox + -"kx + ""L + '2vOa:Vkx) 



a(l -|- eo cos ^o) 



G(A4o - Am) 



which are consistent with iKalogeral (|l996l ). In i ]2.3l we 
present a few examples regarding the effect of mass loss and 
the supernova kick on the orbital parameters of hierarchical 
triples. To compute the systemic velocity of the binary sys- 
tem due to the SN, we begin by writing the pre-SN velocities 
of the primary and secondary in the cm reference frame; us- 
ing the pre-SN mass ratio /io = m2/Mo, these velocities are 
given by 



= fJ'o(^VOx,VOy,()j , 

{fiO - l)(^«0a:,l'0!/,0y 



Vl 

V2 = 



(9) 
(10) 



As a result of the assumption of an instantaneous SN and 
neglecting the shell impact, the instantaneous velocity of the 
secondary remains unchanged after the SN, but the instan- 
taneous velocity of the primary changes to 



^.= {>.OVOx+Vk.,^OVOy+Vky,V,.). 



(11) 



We now use the post-SN mass ratio /i — m2/M, and find 
the systemic velocity of the binary system: 

Vsys = (1 - m)vi + MV2 

= (1 — ^)1-; VOx+Vkx,- VOy +Vky,Vkz 



1-11 



(12) 



These results are consistent with the previously mentioned 
studies on SN in binaries. As a conseqence a binary in which 
the compact object does not receive a kick in the supernova 
explosion moves through space like a frisbee. 

2.1.1 Dissociating binary systems 

The mass loss and the kick velocity have a potentially dis- 
rupting effect on the binary system. However, in cases where 
the mass loss alone would have been large enough to unbind 
the binary, the combination of the two can result in the 
binary system surviving the SN (Hills 1983). If the binary 
system dissociates, the two stars move away from each other 
on a hyperbolic or, in a limiting case, a parabolic trajectory. 
This corresponds to the cases where o < and e > 1 (hy- 
perbola) or a — > oo and e — 1 (parabola). From equation ((Tjl 
we see that for a dissociating binary the angle (p between the 
kick velocit y Vk and the pre-SN relative velocity Vo satisfies 
(|Hillslll983l ); 



1 



2ao Am 
~ MT 



via I \ Vc,o 



2ao 
r 



.(13) 



If the right-hand side of equation (I13p is less than —1, the 
binary dissociates for all i\>\ but if it is greater than 1 the 
binary survives for all 0. If the right-hand side is within the 
range —1 t o 1, the probability of dissociating the binary is 
l|Hillslll983l ): 



2ao Am 
r Mo 



(8b) 



(14) 

iTauris fc Takend (|l998l ) presented analytical formulas to cal- 
culate the dissociation velocities for a binary with a pre- 
SN circular orbit. We follow Tauris & Takens' calculations, 
though ignore the SN shell impact, to derive the runaway 
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velocities of two stars in dissociating binaries, however we do 
so for a pre-SN orbit with arbitrary eccentricity. We use the 
cm coordinate system, explained above. Using the following 
shorthand relations 
M 

Mo 



J = 



2 

Vox 



2m- 



ao 



2ao 



+ 



Vf. 2vOxVkx 



+ 



= 1 + 



j 2ao - r v^y 2ao - r 



m ao 



ao 



1/77 2ao - r j 2ao~r 

- Vky Z 

/i \ mvo ao m ao 



Vky 

fi\ rh ao 



(fc + 1) 



we find the runaway velocities for the primary and secondary 
star: 



(vkx (y + l) + (y + fJ-O^VOx, fJ.OVOy 

tJ,Q^VOx, {^^0 — l)vOy 

(16) 



Vl.dii 



V2,diss = : - : + 1 - 

m2t Vm2t 



,Vky fcVj Ufc^ \ 

H Wo, J j 

m2n m2n m2t/ 



2.2 Hierarchical triple systems 

We now consider a hierarchical system of three stars with 
the primary, secondary and tertiary star having mass, po- 
sition and velocity given by (?ni_o,ri,vi,o), (m2,r2,V2) and 
(m3,r3,V3) respectively. The primary star undergoes a SN 
and the inner binary configuration and parameters are the 
same as in section 12.11 The effective mass of the inner bi- 
nary's centre of mass (cm) is rricmfi ~ rnifi + m2 — Mo, and 
is at position 



rcm,o = (1 — A*o)ri + fJ-0T2 
and has a velocity 

Vcm.O = (1 — Aio)vi,o + AtOV2. 



(17) 



(18) 



The cm and tertiary constitute an outer binary defined by 
the semi-major axis, ^o, eccentricity, Eo, and true anomaly, 
Qo- The separation distance between the cm and the tertiary 
star we denote by Ro. Before the SN the outer binary orbital 
plane has an inclination io with respect to the inner binary 
and the separation distance of the outer binary projected 
onto the xy-plane makes an angle ao with the separation 
distance of the inner binary. This inner-outer binary con- 
figuration is to some extent acceptable, because the triple 
is hierarchical. This implies that the separation distance of 
the cm and the tertiary is large compared to the separation 
distance of the primary and secondary, i.e. _Ro ^ ro, so that 
the tertiary experiences gravitational infiuence of the inner 
binary as if it was comingfrom one star at the cm. We as- 
sume an instantaneous SIs|j. Due to the primary undergoing 

^ See section 12.11 and note that the statements about the inner 
companion (the secondary) also hold for the outer companion (the 
tertiary). 



a SN, the inner binary experiences a mass loss Am and an 
effective kick velocity is imparted to the cm: the systemic 
velocity of the inner binary Vsys given by equation 1121 In 
addition, because of the reduction in mass of the primary, 
the position of the cm has changed due to an instantaneous 
translation along the y-axis 



AR 



(m - Mo) 



ao(l - Co) 



-(o,i,o). 



(19) 



1 -I- eo cos 9o 

The orbital parameters change as a result of the SN: the 
inner binary parameters change according to the description 
in section [2m and the outer binary orbital parameters change 
to semi-major axis, A, eccentricity, E, and true anomaly, 
O. The hierarchical triple before the SN has a total mass 
A'It,o = Mo + ma. We use the cm coordinate system to pin 
down the inner binary and add to this coordinate system 
the tertiary at a position such that Ro ^ ro (see Figure [2]). 
We now select a reference frame in which the center of mass 
of the triple (CM) is at rest (the CM reference frame) . 

Prior to the SN the separation distance between the cm 
and the tertiary is 

Ao{l~Ei) 
1 + -Bo cos 6o ' 
and, using the following shorthand notation 



Ro 



cos io sin ao , — cos io cos ao , sin io ) , (20) 



X 

Y 
X' 
Y' 



Ao\ 1 — Eg cos To cos So + Ao sin To sin 6o 



— ^0 y 1 — -E-o cos To sin Oo + Ao sin To cos Oq 
X cos ao — Y cos io sin ao 
X sin ao +Y cos io cos ao 



Z' = Y sin io 

Vox = Vo 

Voy = Vo 

Vo. = Vo 



X' 










Y' 












z' 






+ Y'2 


+ ^'2 



in which To is the pre-SN outer orbit eccentric anomaly de- 
fined by -Ro ~ ^o(l — -BqcosFo), the velocity of the cm 
relative to the tertiary is 



Vo = Vcm,0 - V3 = {Vox, Voy,Voz)- 



(21) 



The effective kick velocity Vsys makes an angle $ with the 
pre-SN relative velocity of the cm with respect to the tertiary 
star Vo. After the SN the separation distance between the 
cm and the tertiary star is 



R = Ro -f AR, 

Ao(l-E? 



1 -I- -Bo cos 



oL( 

Oo V 



COS io sin ao, {fi — fio) 



ao(l - eg) 
1 -I- eo cos do 



1 -I- -Bo cos So . . . \ 

X -jT cosiocosao,smio , (22) 

Ao(l--Bo) / 

the velocity of the cm relative to the tertiary star is 

V = Vo + Vsys 

~ {Vox + Vsys,x, Voy + Vsys,y,Voz + Vsys, z), (23) 

the cm mass is rricm ~ Mo — Am and the total triple mass is 
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a. The cm coordinate system in the CM reference frame for a 
hierarchical triple system before the SN (at t = 0). 
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b. The cm coordinate system in the CM reference frame for a 
hierarchical triple system after the SN. 

Figure 2. Schematic representation of a hierarchical triple star 
system in the pre- and post-SN phase. The solid blue circles de- 
note the primary and secondary (inner binary); the solid red cirles 
denote the cm and the tertiary (outer binary); the green cirle de- 
notes the CM. The solid arrows denote the velocities the stars or 
cm have at that phase; the dashed arrows denote the velocity the 
SN imposes on the stars or cm which will change its velocity in the 
next phase, (a) In the pre-SN phase at the moment immediately 
preceding the SN, the coordinate system is centered on cm and 
CM is at rest, (b) In the post-SN phase the coordinate system 
is no longer centered on the cm - the cm has been translated in 
the y-direction, towards the secondary - and the CM is no longer 
at rest. In both cases the inner binary orbital plane lies in the 
xy-plane and the y-axis is the line connecting the primary and 
the secondary. 



Mt = Mt.o — Am. The inclination of the outer binary orbital 
plane with respect to the inner binary orbital plane is given 
by: 



IrI 



■ smiQ. 



(24) 



The angle of the outer binary separation distance projected 
onto the xz-plane relative to the inner binary separation 
distance is given by: 

IRol cosio 



IRi 



■ smao- 



(25) 



Applying the relevant equations above and equations ([T]) 
and ((2)1 to our triple system, we obtain equations relating the 



post-SN semi-major axis. A, and eccentricity, E, to both the 
pre- and post-SN orbital parameters and velocities. Using 
Vc,o = V'olflo^Ao = {GMt,o/Ao)^^^ as the pre-SN relative 
velocity when Ro = Ao, and using p = (_Ro — R)/{RoR), we 
obtain 



Ao V A/t,o/V 



Am - 



2Ao Am 



2— — — ^ cos$ 



R Mt,o Vc,o Vc,o 



+ 2Aop 



(26) 



1 - (1 - El] 



Mt,o 



(2Ao 



Mt.o 



xl 



2Ao 
Ro 



{Mt,o - Am) V R 

«.2 



Mt,o - Am 



-»M» r, Vsys -,\\ 

Tjt - 2t^T7^cos$ . 



(27) 



With the pre-SN mass ratio vo ~ mz/Mt,o, the pre-SN ve- 
locities of the cm and the tertiary in the CM reference frame 
are 



Vcm.O = Vo(Vox,Voy,Voz^ 

(uo - 1) (Vox,Voy, Vo, 



V3 = 



(28) 
(29) 



We calculate the instantaneous velocity of the cm after the 
SN (as before, because of the assumption of an instanta- 
neous SN, the velocity of the tertiary after the SN remains 
unchanged) : 



= fO Vox + 



'Vsys,x ,r Vsys,y Vsys,z 
Voy H , VOz + 



Vo 



Vo 



(30) 



Using the post-SN mass ratio u = mz/Mt, the systemic 
velocity of the outer binary (and therefore of the triple) is 



(1 - i^)Vcm + vwa 

Vo - V 



vo — V 
1 - V 

vo — V 
1 - V 



Vox + (mo — ^J^VOx + (1 — p)vkx, 



1 - V 

Voy + (mo — H)'VOy + (1 — \i)Vhy, 
Voz + (1 - H)Vk^- 



(31) 



Summarizing, one can consider a hierarchical triple system 
as a effective binary system composed of an effective star 
(i.e. the inner binary center of mass (cm)) and the tertiary. 
The effective star undergoes an effective asymmetric SN re- 
sulting in three effects: 1) sudden mass loss Am, 2) an in- 
stantaneous translation AR, and 3) a random kick velocity 
Vsys. The calculation of the post-SN parameters and veloc- 
ities of a hierarchical triple system is now reduced to the 
prescription for a SN in a binary as presented in section [2TT] 
Note that the mass loss does not occur from the position of 
the effective star, but from the position of the primary star; 
a clear distinction from a physical binary system. However, 
from what position the mass loss occurs is not important 
when an instantaneous SN is considered. When the effect 
of the shell impact on the companion star(s) is considered, 
this off-center mass loss must be taken into account. In ad- 
dition, if it were not the primary which underwent the SN, 
but, for example, the tertiary, the computation would have 
been done by reducing the inner binary to an effective star, 
as shown in this section. One would again have a binary con- 
figuration to calculate the effect of the SN; in such a system 
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there is no off-center mass loss. In section \2A\ we show how 
one can reduce any hierarchical multiple star system to an 
effective binary in a recursive way using the effective binary 
method and in § 12.4.31 we do the computation of the effect 
of a SN on a binary-binary system. 



2.2.1 Dissociating hterarchical triple systems 

For the triple system, dissociation can occur in two ways: 
the inner binary can dissociate (a < and e > 1 or a — >■ oo 
and e = 1) (see section [27111 and the outer binary can disso- 
ciate < and iJ > 1 or A — ^ oo and E = 1), i.e. the inner 
binary and the tertiary become unbound. The inner binary 
dissociation scenario generally results in complete dissocia- 
tion of the system. However, hypothetical scenarios exist in 
which one of the inner binary components is ejected towards 
the tertiary star to either collapse with it or to form a bi- 
nary by gravitational or tidal capture. Nevertheless, these 
scenarios have a small probability since the ejection con- 
ditions (e.g. the solid angle in which that particular inner 
binary component has to be ejected in) and the capture con- 
ditions are extremely specific. From equation[26]we see that 
for the inner binary to dissociate from the tertiary, the angle 
$ has to satisfy 



cos $ > 



2Ao Am 



\/2 



2A,p)(2 "^'^ 



Vc,0 Vc,0 



(32) 



The probability of this type of dissociation is 

pouter 



1 ( ( 2ylo Am vl \ 
Vo V. 



X 2 



Vc 



(33) 



In the case of the dissociation of the outer binary, using the 
following short hand relations 



M 



Mt 
Mt.< 



K 



^ ^ J 2Ao-Ro R _ vlys,y 2Ao-R R 
M Ao Ro MV(f Ro 

1 f VJ 2Ao- Ro R 

-Vs 



i^^MVo""^"'^ ^0 Ro 
J 2Ao -Ro R \ 
M Ao Ro ) 

N = lfi + ^2^°-^°^(j^ + l)) 
v\ M Ao Ro I 

the runaway velocities of the inner binar y system and the 
tertiary are (following and generalizing iTauris fc Taken^ 



|T99i)): 



Vcm.diss = \v sy s, x\ + + VO^Vox.Vsy s,y{^— 

+ T^oVoy + Vo,Vsys,z(^J- + l) ) (3'? 



V3 Hiss — 



N 
/ 1 



m-iL \m:jL 



+ 1 — uoj Vox, iyo — '^)Voy 



Vo, 



m:jN rri'iN ' msL 



(3E 



Note that these equations are more general than the ones 
in section r2.1.1l because we cannot assume R = Ro in the 
triple case. 



2.3 An example of the effect of a supernova in a 
hierarchical triple 

For two simple sets of initial conditions we investigated the 
effect of mass loss. Am, and kick velocity, Vk, on the surviv- 
ability of a triple system. We distinguish between four dif- 
ferent post-SN scenarios: (1) the triple survives as a whole 
(e < 1 and _B < 1) with new orbital parameters, (2) the 
inner binary survives and the third star escapes (e < 1 and 
E > 1), (3) the inner binary dissociates and the outer binary 
survives (e > 1 and E < 1) and (4) the triple completely 
dissociates (e > 1 and E > 1). The third scenario is a rather 
special case and can only be of temporary nature: in this 
scenario, even though the inner binary has just dissociated, 
the third star remains bound to the inner binary center of 
mass. This is a temporal solution which eventually will lead 
to the full dissociation of the triple, except in the extreme 
case in which the tertiary star captures one of the ejected 
inner stars to form a new binary system. 

For each set of initial conditions we used a hierarchi- 
cal triple system with primary, secondary and tertiary stars 
of masses mi^o, ni2, ma = 3, 2, 1 Mq respectively and in- 
ner and outer binary semi-major axes ao, Ao = 10, 50 Rq 
respectively, and we varied the kick velocity direction Vk. 
For the two different sets of initial conditions we determine 
which combinations of Am and Vk lead to which post-SN 
scenario and we show our results in Figure O the used ini- 
tial conditions are specified below the respective figures. 

In Figure [3^. we used a circular inner and outer orbit, 
not inclined with respect to each other, with all stars on 
one line and the kick velocity in the same direction as the 
pre-SN inner binary relative velocity. We see that for zero 
kick velocity, the inner binary dissociates for a mass loss 
ratio of Am /Mo = 0.5, which is consistent with earlier work 
(e.g. iHill j [l983l ) . For zero mass loss, we see that the inner 
binary dissociates for a kick velocity of Vk ~ 128 km/s - this 
velocity is exactly the difference between the inner binary 
escape velocity {vesc = ^j2GMol ao 
SN relative velocity (vo 
the third star escapes for a slightly lower value of the kick 
velocity. This is because the inner binary systemic velocity 
(which is the effective outer orbit kick; see Section [2. 2p plus 
the pre-SN outer orbit relative velocity already exceed the 
outer orbit escape velocity. We furthermore see that the total 
triple survival scenario allows lower kick velocities for higher 
mass losses. Above a kick velocity of Vk ~ 128 km/s the 



437 km/s) and pre- 
^jGMo/ao ~ 309 km/s) - but 
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inner binary always dissociates, irrespective of the mass loss, 
(eventually) leading to total dissociation. 

In Figure [SJj- we keep the same configuration as de- 
scribed for Figure |3^., but with a kick velocity in the op- 
posite direction with respect to the orbital velocity of the 
exploding star before the supernova. The triple can now lose 
more mass and receive a higher velocity kick while stil sur- 
viving. The ability to sustain greater kick velocities is ex- 
plained by the fact that, depending on the mass loss, the 
kick velocity now has to exceed a fraction of the sum of vo 
and Vfc (for zero mass loss vo-|-Vfc~ 746 km/s) due to the 
opposing directions of the two velocities. We also see that 
total triple survival can occur beyond a mass loss ratio of 
0.5, because the kick velocity can opp ose the dis sociating 
effect of the mass loss (as mentioned in lHills|[l983l ). Bear in 
mind that while the Am /Mo — case is non-physical we 
include it for the sake of completeness. 

In Figure[3]we show how the post-SN systemic velocity 
of the triple depends on the mass loss Am for a hierarchi- 
cal triple system with primary, secondary and tertiary stars 
with masses (mi_o, m2, ma) — (3, 2, 1) Mq, inner and outer 
binary semi-major axes (ao, Aq) = (10, 50) Kq and the kick 
velocity in the direction of the pre-SN inner orbit relative 
velocity. We plot our results for the case that the SN went 
off at the inner orbit apastron (6o — 180 degrees) or at the 
inner orbit periastron {9q — degrees) for a symmetric SN 
(i.e. Vk = km/s) and a SN with a kick Vk ~ 31 km/s, in the 
cm reference frame (i.e. with the cm at rest at t = 0). In the 
top panel of Figure[4]we see that for a symmetric supernova, 
the systemic velocity of the inner binary increases with the 
amount of mass loss, which is an intuitive result. We see that 
even with zero mass loss the triple has a systemic velocity, 
namely the velocity it started with in this reference frame 
{Vays ~ 17.5 km/s). We furthermore see that the increase of 
the triple systemic velocity happens more steeply for these 
cases where the SN goes off at periastron - with the steepest 
curve for the highest inner binary eccentricity - than when 
the supernova goes off at apastron - with the steepest curve 
is for lowest eccentricity. For an asymmetric supernova with 
kick Vk ~ 31 km/s (see the bottom panel of Figure[4]) we ob- 
serve similar behaviour, but with the difference of the zero 
mass loss case: in this case the triple 

system has a lower velocity than it started with {Vsys ~ 
2.5 km/s), which is due to the kick. This result is dependent 
on the direction of the kick. 

The pre-SN triple systemic velocity is dependent on both the 
inner binary and the outer binary. Its dependence on the in- 
ner binary is via the masses mi.o and m2 of the primary 
and secondary respectively and the inner binary orbital pa- 
rameters which fully constrain the relative velocity of these 
stars (see equation (O). Its dependence on the outer binary 
is via the mass ms of the tertiary and the outer orbit orbital 
parameters which fully constrain the outer binary relative 
velocity (see equation (|2ip l. The post-SN triple systemic 
velocity is merely the sum of the pre-SN systemic velocity 
and its change, which is only due to the inner binary through 
the mass loss Am and kick velocity Vk. 

2.4 Hierarchical systems of multiplicity > 3 

There exist two kind of hierarchical multiple star systems 
with more than three stars: 




%.0 0.1 0.2 0.3 0.4 0.5 

Am/Mg 




!l.O 0.1 0.2 0.3 0.4 0.5 



Am./M„ 

Figure 4. The post-SN systemic velocity of the triple as a func- 
tion of mass loss. Am, when the SNe occurs at periastron {Bq = 0, 
dashes) and apastron (9o = 180 degrees, solid curves) of the inner 
binary, for a range of pre-SN inner binary eccentricities, = 
km/s in the top panel and ^ 31 km/s in the bottom panel. 

(i) systems that have n stars and hierarchy n — 1, i.e. 
multiple star systems with its stars hierarchically ordered in 
series (hereafter serial systems). Examples of such systems 
include quadruples with hierarchy 3, but also binaries and 
triples are serial systems. 

(ii) systems that have n stars and hierarchy n—2 or below, 
i.e. multiples composed of serial systems which are hierar- 
chically ordered in parallel (hereafter parallel systems). An 
example of such system is a quadruple with hierarchy 2 (i.e. 
a binary- binary system). 

2.4-1 Serial systems 

The effect of a SN on a serial system is calculated by ap- 
plying the effective binary method (see section 12. 2p by re- 
cursively replacing the inner binary by an effective star at 
the center of mass of that binary, until the total system 
is reduced to a single effective binary. When considering a 
serial system of n stars each with mass, position and ve- 
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A)n,/Af„ Am/Mf. 



(a) eo = 0, Eo = 0, do = 0°, 0o = 0°, io = 0, «o = 0°, (b) eo = 0, Eq = 0, 6*0 = 0°, 60 = 0°, io = 0, «o = 0°, 

Vk = (1,0,0) Vk = (-1,0,0) 

Figure 3. The plots above show the survivabihty of the hierarchical triple system for varying mass loss Am and kick velocity v^- The 
systems have masses of mi^o, ni2, ma = 3, 2, 1 Mq respectively and inner and outer binary semi-major axes ao, Aq = 10, 50 R© 
respectively. There are four possible post-SN scenarios: (1) the whole triple survives, (2) the inner binary survives but the third star 
escapes, (3) the inner binary dissociates and the outer binary survives, or (4) the triple completely dissociates. The areas in the plots 
are labeled according to their respective post-SN scenario. 



locity given by (mi,o,ri,Vi,o), (m2,r2,V2), ... , (m„,r„,Vn) 
respectively, in which the primary star undergoes a SN, one 
starts by reducing the inner binary to an effective star, as 
was done in section 12.21 The inner binary consists of the 
primary and secondary star at positions ri and r2 respec- 
tively. This binary is reduced to an effective star of mass 
TTicmfi = JTii^o + ?Ti2 at position Tcm.o giveu by equation l|17p 
and having velocity Vcm.o given by equation (|18|l . Due to 
the SN of the primary this effective star experiences a mass 
loss Am, an instantaneous translation AR given by equa- 
tion (|19p , and a random kick velocity Vsys given by equation 
(|12|) . After applying these effects on this effective binary, one 
can calculate the post-SN orbital parameters and velocities 
and the systemic velocity Vsyi = Vsys of this effective bi- 
nary, given by equation (|3ip . using the prescription for a SN 
in a binarylj The total system is now reduced to a serial 
system of n — 1 objects (real and effective stars) . 

Subsequently, one reduces the current inner binary - 
consisting of the effective and tertiary star at positions rcm.o 
and ra respectively - to an effective star of mass m).^ q = 
m-cm,o + 1^3, at position 



^(2) _ mcn.prcm.O + mzTs 

with a velocity 

(2) T^cm, oVcm.O + maVa 

rricmfi + ma 



(36) 



(37) 



Due to the SN of the primary star, this effective star also 
experiences a mass loss Am, an instantaneous translation 
AR^^^ - this time, the translation vector has non-zero y- 



^ The number between parentheses denotes the hierarchy up to 
which the system has been reduced to a effective star. 



and z-components - and a random kick velocity Vsys. Af- 
ter applying these effects on this effective binary, one can 

calculate the post-SN orbital parameters and velocities and 
(3) 

the systemic velocity Vsys of this effective binary using the 
prescription for a SN in a binary. The total system is now 
reduced to a serial system of n — 2 objects (real and effective 
stars). 

This procedure is carried on until the entire multiple 
is reduced to a single effective binary, consisting of the nth 



star at position rn and a effective star of mass m^ 
m^^~o' -|-m„_i at position 



(n-2) 



' cm,0 



(n-3) {n-3) 
(1-3) 



with a velocity 



,(n-2) 



(n-3) (n-3) 
'cm,0 "crn.O 



+ m„_lVn-l 



(n-3) 



(38) 



(39) 



This effective star also experiences mass loss Am, an instan- 
taneous translation AR'-""'^^ and a random kick velocity 
viy^'^'- After applying these effects on this (final) effective 
binary, one can calculate the post-SN orbital parameters and 
velocities and the systemic velocity viy^"""' for this effective 
binary (and therefore of the total system) using the binary 
method. 

When it is not the primary star which undergoes a SN, 
but the mth star in the hierarchy, the procedure is carried 
out by first reducing the inner serial system of m — 1 stars 
to an effective star at its center of mass. One can then apply 
the above explained method, as there is no computational 
difference in whether the primary or the secondary of a(n 
effective) binary undergoes the SN. 
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2.4-2 Parallel systems 

The effect of a SN on a parallel system is calculated by 
reducing each parallel branch (which itself is a serial system) 
to an effective star until an effective serial configuration is 
reached; after this, one can use the method explained in the 
previous section. We consider a parallel system of i parallel 
branches, each consisting of an arbitrary number Ui of stars 
with mass, position and velocity given by (mi,ri,vi), ... , 
(m„; ,rni ,Vni) respectively, in which the mth star - which is 
part of branch j - undergoes a SN. One starts by reducing 
alH — 1 branches 7^ j to effective stars. One then calculates 
the effect of the SN on branch j (i.e. systemic velocity and 
mass loss) using the method described in section [2.4.1l The 
total system is now reduced to an effective serial system of 
i effective stars in which the jth effective star undergoes 
an effective SN with the systemic velocity of branch j as 
the kick velocity. The effect of this effective SN on the total 
system, can be calculated by applying the method described 
in section [2.4.1l to this effective serial system. As an example 
we will now demonstrate the effect of a SN on a binary- 
binary system. 

2.4-3 An example of the effect of a supernova in 
binary-binary system 

We consider a hierarchical binary-binary system of stars 
with mass, position and velocity given by (mi,o,i'i,vi,o), 
(m2,r2,V2), (m3,r3,V3) and (m4,r4,V4) respectively, in 
which the primary star undergoes a SN. The binary consist- 
ing of the primary and the secondary star (primary binary) 
has the configuration and the parameters as in section 12.11 
and has a center of mass (cmi, i.e. effective star 1) of mass 
m-cmi,o = JTii.o + m2 = Mo at position given by equation 
(|17p with a velocity Vcmi,o given by equation p8|) . The sec- 
ondary binary consists of the tertiary and quaternary star 
and its center of mass (cm2, i.e. effective star 2) has a mass 
mcm2 ~ 'TIS + rn4 = M2, is at position 

rcm2 = (1 - K)r3 + Kr4 

and has velocity 

Vcni2 = (1 - «:)V3 + KV4, 

before the SN, where k = The cmi and cm2 constitute 
an effective binary defined by semi-major axis, Aq, eccen- 
tricity, Eo, and true anomaly, Oo. The separation distance 
is denoted by Rq. Before the SN the effective binary orbital 
plane has inclination io with respect to the primary binary 
orbital plane and the separation distance of the effective bi- 
nary projected onto the xy-plane makes an angle ao with the 
separation distance of the primary binary. We assume an in- 
stantaneous In the effective SN the cmi experiences a 
mass loss Am, an instantaneous translation AR along the 
X-axis given by equation (|19|) and a random kick velocity 
Vsys given by equation p2|l . The orbital parameters change 
as a result of the SN: the primary binary parameters change 
according to the description in section 12.11 and the effec- 
tive binary orbital parameters change to semi-major axis A, 

* See section [2. II and note that these statements about the inner 
companion (secondary) star also hold for the outer companion 
(tertiary and quaternary) stars. 



eccentricity E and true anomaly O; the secondary binary or- 
bital parameters do not change when SN-shell impact is not 
taken into account. Before the SN the binary-binary system 
has a total mass Mbb,o = mcmifi + mcm2, we use the cmi 
coordinate system to pin down the primary binary and add 
to this coordinate system the tertiary and quaternary at a 
position such that 7?o S> ^o, and we choose a reference frame 
in which the center of mass of the total binary-binary system 
(CMbb) is at rest (the CM;,;, reference frame) and in which 
the cmi is at the origin at t — 0. The separation distance 
between the cmi and the cm2, Ro, is given by equation (|20|l 
and the velocity of the cmi relative to the cm2 is 

Vo = Vcml.O - Vcm2 = {Vox,Voy, Vqz) (40) 

prior to the SN. The effective kick velocity Vsys makes an 
angle $ with the pre-SN relative velocity Vo. After the SN 
the separation distance between the cmi and the cm2 is R 
given by equation (|22|l and the velocity of the cmi rela- 
tive to the cm2 is V given by equation (|23|l . the cmi mass 
rricrni = rncmi,o — Am = M and total binary-binary mass 
Mtb = rricmi + rncm2 — M + M2. Applying the relations 
above and equations ([1} and ([2]) to our binary-binary sys- 
tem, we obtain relations for the post-SN semi-major axis A 
and eccentricity E in terms of both the pre- and post-SN or- 
bital parameters and velocities given by equations (|26p and 
(|27|l respectively with Mtfi replaced by Mbb,o- To compute 
the systemic velocity due to the SN, we express the pre- 
SN velocities of the cmi and the cm2 in the CMj,;, reference 

frame. Using the pre-SN mass ratio An ~ tt^, the pre-SN 

° -''-'66,0 

velocities are given by 

Vcml.O = \Q(Vox,Voy,VQ2^ (41) 

Vcm2 = (Ao-l)(M),,Voy,l/o,). (42) 

We calculate the instantaneous velocity of the cmi after the 
SN (due to the assumption of an instantaneous SN, the ve- 
locity of the cm2 after the SN remains unchanged): 

vc^i = Ao (Vo. + ^^,Voy + Vo. + ^) (43) 

V A() Ao Ao / 

With the post-SN mass ratio A — , the systemic veloc- 

ity of the effective binary (and therefore of the binary-binary 
system) is 

Vsys = (1 - A)Vcml + AVcm2 

= (1 — A) ^ Y — ^^'^ + ^ tJ-)vox + (1 — ^J.)vkx, 

y— ^Voy + (mo - I^)y0y + (1 - li)Vky, 
^^^Vo, + {l-IJ.)vk.)- (44) 

Note that because the branch harboring the SN-progenitor 
(SN branch) is a binary, this calculation the SN-effect on 
the binary-binary system is almost identical to calculation 
of the SN-effect on a hierarchical triple. The computations 
become more interesting for systems with a SN branch of 
higher multiplicity. 
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3 APPLICATION: FORMATION OF 
J1903+0327 

PSR J1903+0327 was observed by IChampion et all (|2008l l 
who determined it to be a millisecond pulsar (MSP). This 
MSP is observed to have a 1 Mq main sequence compan- 
ion with a highly eccentric and distant orbit (e ~ 0.44, or- 
bital period ~ 95.2 days). These properties are atypical for 
MSPs be cause MSPs are expected to be spun-u p via mass 
transfer (jBhattacharva fc van den Heuvell[l99ll ). which in 
turn widens and circularizes t he orbit, while its companion 
evolves through a giant phase. iPhinne \1 (|1992| '). for example, 
suggest an eccentricity e < 10~^ is typical for MSP binaries. 
The exception to this has been MSPs in globular clusters 
which have interactions with other obj ects that m a y per - 
turb the orbit of the binary. However, iFreire et al.1 (l201ll ) 
find it to be unlikely that this MSP system has its origin in 
an exchange interaction in such a dense stellar environment. 

It has been suggested that J 1903-1-0327 may b e the 
result of a hierarchical tripl e ([Champion et al.l |2008| . 
iPortegies Zwart et al.l l201ll and iBeiger et al.1 120111 ) where 
the inner companion has been lost after spinning-up the 
MSP, leaving only the MSP and the former tertiary to be 
observed. Should J1903-I-0327 be the result of such a sys- 
tem the methods in the previous sections provide a strong 
beginning to investigate how such a system might evolve. 



3.1 Initial conditions 

We generate sets of 10^ initial conditions, as described be- 
low, with each set constituting a stable triple system, and 
then simulated the effect of an instantaneous SN occurring 
at the primary star. The mo del we follow (many of our in i- 
tial conditions are drawn from lPortegies Zwart et ahl (|201ll )) 
consist of a primary, secondary and tertiary star with zero 
age masses of 10 Mq, 1 Mq and 0.9 Mq respectively. The 
initial conditions are generated by selecting the semi-major 
axis, Ao, eccentricity, Eo, and the orbital inclination, i, for 
the tertiary. Aq takes values on the range [200, 10 000]i?Q 
from a flat distribution, Eq is chosen on the range [0, 1) from 
a distribution that is flat in log space, and io is chosen on 
the range [0, tt] with a sinusoidal distribution. Combining 
these values with the zero age masses of the stars as well as 
a pre-set value for the initial semi-major axis of the inner 
binary, ao = 200i?Q we then test for stability of the system 
using: 



> 3 14--— - + - cos Jo - cos 10 

oo V Mq/ V4 2 



1/3 



En 



-1/6 



(45) 



l|Zhuchkov et al.ll2010l ). If the system is stable with this set 
of parameters, we choose the remaining parameters, namely 
the angle ao described in the previous sections, the direction 
and magnitude of the kick. Because we have assured that the 
system is dynamically stable before starting our simulations 
our assumption of a hierarchical system is guaranteed. We 
observe that due to the SN kick, systems with very high 
inclination are preferentially removed or their inclination is 
reduced thus as a result we do not include the effects of 
Kozai iterations. 




60 80 
Vk (Km/s) 

Figure 5. The fraction of surviving and stable system (thin red 
and thick black (colors online only) lines respectively) as a func- 
tion of the kick velocity. The lines in each set correspond to differ- 
ent semi-major axis, 50, 30, 20, and 10 Rq (circle, cross, diamond, 
and square respectively). All curves are normalized to the total 
number of surviving systems with a semi-major axis of 50 Rq. 



3.2 Simulations 

The inner binary undergoes a common envelope (CE) phase, 
circularizing the orbit, reducing the inner semi-major axis to 
a value between 5 Rq and 60 Rq , and reducing the mass of 
the primary to 2.7 Mq. The effect of these changes on the 
stability of the system can immediately be seen in equation 
(|45|l . Then, due to the SN, the primary undergoes a mass 
loss of 1.3 Mq and receives a corresponding kick. The ve- 
locity of the kick is fixed between 5 and 160 km/s for each 
set of simulations and the kick direction is randomly chosen 
such that for all simulations the direction is isotropic. We 
then analyze the survivability and stability of each system. 
A system survives the SN and resulting kick if it remains 
bound, and it is determined to be stable if, while remain- 
ing bound, the system also satisfies the stability criterion in 
equation (|45|l . 

We ran Monte Carlo simulations for four different inner 
binary semi-major axes (10, 20, 30, and 50 Rq). For each 
semi-major axis value we run 25 simulations (each of the 
25 simulations consists of 10^ sets of initial conditions) each 
with a constraint kick velocity (between and 130 km/s). 
In Figure Owe plot the kick velocity versus the fraction of 
surviving and stable systems. For each pair of curves the 
thin red upper curve corresponds to the survivability frac- 
tion and the thick black lower curve to the fraction that 
survives and remains stable. Curves with same kick veloc- 
ity have the same point-symbols. Each point represents the 
fraction of surviving or stable systems normalized to the 
total number of surviving systems with a semi-major axis 
of 50Rq. Increasing the semi-major axis from 10 to 30 Rq 
strongly increases the overall probability of a system to sur- 
vive and remain stable. However, with a kick velocity of 
45 km/s and higher the probability of a system remaining 
stable is nearly the same when the semi-major axis is ^ 
20Rq . Figure [5 ] shows the effect of the Blaauw & Boersma 
iiwl 1 



recoU (|BlaauwE96il &|b oersmal (|l96ll )) on the system when 



the SN kick is small; as the SN kick velocity approaches the 
Blaauw & Boersma recoil velocity the stability increases due 
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Figure 6. The fraction of surviving and stable systems (upper 
and lower lines respectively) with respect to the inner semi-major 
axis. A constant kick velocity of 20 km/s is used. 
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Figure 7. The number of occurrences for which the system be- 
comes unstable due to mass transfer at a given mass of the pri- 
mary. The curves corresponds to Face values of 0.3, 0.4, 0.5, 0.6, 
and 0.9 as shown in the key. The peak value and FWHM for each 
curve in this figure, as well as similar curves for other values of 
Face, are plotted in Figure |8] 



to the kick and recoil off-setting one another, in part or in 
full. As the SN kick velocity increases it begins to overwhelm 
the Blaauw & Boersma effect. 

In Figure [6] we show the effect the inner semi- major 
axis has on survivability and stability (the upper and lower 
lines respectively) using a constant kick velocity of 20 km/s. 
Again each data point represents the fraction of systems that 
survive or survive and in addition remains stable out of a set 
of 10^ initial conditions. Hfere we see the significant role of 
the inner semi-major axis on the survivability of the system. 
If we note for a particular kick velocity which value of ao the 
stability fraction begins to level, we can see it corresponds 
to the merging of the stability curves in Figure [S] For the 
case of a 20 km/s SN kick velocity, as in Figure [S] we see 
that any value of ao greater than about 30 Rq will have 
similar stability fractions while systems with lower values of 
ao should have a lower stability fraction as we see in Figure[Sl 

Next, we chose all of the systems that remain stable 
after the SN and subject them to a mass transfer phase. 



Figure 8. The final mass of primaries with respect to the fraction 
of accreted mass. The dashed horizontal line is placed at the ob- 
served mass of J1903-I-0327. The points represent the peak value 
of curves which plot the number of times a system becomes unsta- 
ble while at a given mass of the primary (like those in Figure [7]|; 
the upper and lower bars represent the FWHM of the curves. 
The values that are colored (online) and that have different line 
types correspond to the curves in Figure[7](e.g. the blue, dot-dash 
line at Facc=0.9 is obtained from the right most peaked curve in 
Figure [T] which is also a blue, dot-dash line). 



Here we iteratively remove one one-hundredth of the mass 
of the secondary and transfer a fraction of it to the primary, 
which after the SN would ha ve formed a neu tron star (NS). 
Following the work of Pols fc MarinusI ij 19941 ) we find: 

— [(Sf)'"""""5f]"'KS) 

where a/ is the new semi-major axis, a^ is the semi-major 
axis before the mass transfer, nii^i and m2,i are the masses 
of the primary and secondary before the mass transfer and 
mi J and m2j are the masses of the primary and secondary 
after the mass transfer, Mi and Mf are the total masses of 
the binary before and after the mass transfer, and finally x 
is the ratio of the change in mass of the system to the change 
in mass of the donor (i.e. the secondary). If we define the 
fraction of mass accreted. Face, as the fraction of mass lost 
from the secondary which is accreted onto the primary we 
find that the 1/(1 — x) term simply becomes 1/Facc- After 
each iterative mass transfer, and the resulting change in the 
semi-major axis, we test the triple for stability using equa- 
tion (|45p . When the system becomes dynamically unstable 
we stop simulating as the assumption of a hierarchical sys- 
tem has broken down. We record the mass of the primary 
when the system becomes dynamically unstable and plot 
the mass in Figure [7| versus the number of times systems 
becomes unstable at that mass. For this plot we used Face 
values of 0.3, 0.4, 0.5, 0.6 and 0.9, which correspond to the 
lines which peak from the left to right respectively, and a 
constant kick velocity. We see that the peak value for each 
Face shifts to a larger primary mass as Face increases. This 
relation is expected since as Face becomes larger more of the 
mass lost from the secondary is accreted onto the primary. 
So for the case of Face = 0.3 only 30% of the mass lost from 
the secondary could ever accrete onto the primary thereby 
reducing the maximum possible mass of the primary. If we 
assumed that all of the mass of the secondary is lost (an un- 
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Figure 9. The number of systems per millon simulations with 
a final primary mass of 1.667 Mq (the observed mass of 
J1903+0327) as a function of the fraction of accreted mass, for 
different initial primary masses (shown in the key). 



physical case since the mass transfer would end before this 
could happen, but this provides an extreme upper limit) 
then while the secondary would have lost IMq the primary 
would have only accreted O.SM© resulting in a maximum 
primary mass of 1.7 Mq. If we were to assume that mass 
transfer would stop when the secondary decreased to a mass 
of O.SMq then the secondary would have lost O.7A/0 and 
only 0.21Mq (or 30% of O.JMq) would have been accreted 
by the primary resulting in a mass of 1.61Mq. We have ex- 
amined 21 curves like those in Figure [T] we measured and 
plotted their peak value and the full-width-half-maximum 
(FWHM) in Figure \E\ The error bars denote the FWHM 
of the curves, the plotted point is the peak value for each 
curve, and the mass of J1903-I-0327 is shown as a dashed 
line. Examination of Figure [S] shows that given the observed 
mass and the assumptions we used in preparing the simu- 
lated systems, J1903+0327's progenitor system would have 
most likely had an Face value between between 0.35 and 0.5, 
with the peak value of 0.4 most closly maching the observed 
mass. 

It should be noted however, not all of the barionic 
mass transfered results in an equivalent increase in gravi- 
tational mass of the prima ry since Maccrete ~ AMgrai> + 
^Ehindingl<? (|Bagchill2oTll ) . where Maccrete is the mass ac- 
creted from the secondary, AMgrav is the change in gravi- 
tational mass of the primary, and AEbinding is the binding 



energy of the system. We find that for the masses being 
transferred in our simulations the effect of using Maccrete = 
AMgrav is less than the uncertainty in the final results. 

Finally, we preform the same analysis that produced 
Figure [7] but use an initial primary mass of 1.2, 1.3, 1.4 
(as used in all of the previous simulations), 1.5, and 1.6 
Mq . These simulations were preformed for eight inner semi- 
major axes (10, 20, 30, 40, 50, 60, 70, and 100 Rq) at the 
start of mass transfer. The Face value with the peak number 
of occurrences closest to the observed mass of J1903-I-0327 
(1.667 Mq) was recorded, as was the number of occurrences 
at that peak; these values were plotted in Figure (9] Upon 
examining Figure [5] we find that as the initial mass of the 
primary increases the most likely Face value and its domain 
decrease. To understand these results we recall that as the 



initial mass of the primary increases the amount of mass 
needed to reach the observed mass of J1903-I-0327 is de- 
creased. So, for example, if the initial mass of the progenitor 
of J1903-|-0327's primary (before it began to accrete mate- 
rial from the secondary) was 1.6 Mq it would only need to 
accrete 0.067 Mq before the system reached the observed 
mass. A very small Face value can result in the transfer of 
such a small amount of material allowing the Face to stay 
low; with a lager Face value the system will often reach a 
final primary mass greater than 1.667 Mq thus limiting the 
domain. Whereas if the initial primary mass was 1.2 Mq, 
an Face value of 0.1 would never allow for enough mass to 
be transfered, but there are a large range of Face values that 
can allow for that amount of mass transfer that would not 
quickly overshoot the observed mass. This assumes, as we 
have in all of the simulations, that the mass transfer is stable 
as long as the triple is dynamically stable. We find that for 
an initial primary mass of 1.4 Mq, the value used in all pre- 
vious simulations, the peak Face value is not sensitive to the 
semi-major axis at the beginning of the mass transfer; the 
Face value ranges between 0.35 and 0.45 which lies within 
our expected range of 0.35 to 0.5 found above from Figure|8l 



4 CONCLUSION 

We have examined the effect of an asymmetric supernova 
(SN) on a hierarchical multiple star system and considered 
how it can be modeled by applying the effective binary 
method. This is done by recursively replacing the inner bi- 
nary by an effective star at the center of mass of that binary. 
The effective star experiences an effective SN with the effects 
of sudden mass loss, an instantaneous translation and an ef- 
fective kick velocity, i.e. the systemic velocity of the inner 
binary. We have coded the equations in this paper in a small 
python script, which is publicly availablslf] 

We point out that the effective SN is different from a 
physical SN in that for a physical SN themass is lost from 
the position of the physical star, whereas for an effective SN 
the mass is lost from the effective star. The off-center mass 
loss in an effective SN becomes important only if the shell 
impact on the companion(s) is considered, and otherwise 
causes no difference between a real and effective SN calcula- 
tion. Furthermore, we calculated the runaway velocities for 
dissociating binaries and effective binaries. We subsequently 
demonstrated how calculating the effect of a SN on a mul- 
tiple can be generalized to multiples in which a star other 
than the primary is undergoing the SN. 

We used this method to examine the case for 
J1903-I-0327 forming from a hierarchical triple. We assume 
initial masses of 10, 1.0, and 0.9 Mq for the primary, sec- 
ondary, and tertiary respectively, as well as an inner semi- 
major axis of 200 Rq. We find that if J1903-I-0367 was to 
form through such a mechanism it would be most likely to 
have a very low SN kick velocity so that it would remain 
stable after the SN, and a large inner semi-major axis af- 
ter the CE phase to increase the likelihood that the triple 
would become unstable once the NS/MSP reached a mass of 



^ The source code is publicly available 
|http : //castle ■ strw. leldenuiLlv.nl/sof tware .html | 



at 
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1.667 Mq (|Freire et al.ll201ll '). We also find that, given our 
assumptions, the transfer efliciency, Face, for J 1903+0327 
would have likely been between 0.35 and 0.5. 
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